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ABSTRACT 


The Closeby Habitable Exoplanet Survey (CHES) is dedicated to the astrometric exploration for 
habitable-zone Earth-like planets orbiting solar-type stars in close proximity, achieving unprecedented 
micro-arcsecond precision. Given the elevated precision, thorough consideration of photocenter jitters 
induced by stellar activity becomes imperative. This study endeavors to model the stellar activity 
of solar-type stars, compute astrometric noise, and delineate the detection limits of habitable planets 
within the astrometric domain. Simulations were conducted for identified primary targets of CHES, 
involving the generation of simulated observed data for astrometry and photometry, accounting for 
the impact of stellar activity. Estimation of activity levels in our samples was achieved through 
chromospheric activity indices, revealing that over 9096 of stars exhibited photocenter jitters below 1 
uas. Notably, certain proximate stars, such as o Cen A and B, displayed more discernible noise arising 
from stellar activity. Subsequent tests were performed to evaluate detection performance, unveiling 
that stellar activity tends to have a less pronounced impact on planetary detectability for the majority 
of stars. Approximately 95% of targets demonstrated a detection efficiency exceeding 80%. However, 
for several cold stars, e.g., HD 32450 and HD 21531, with the habitable zones close to the stars, a 
reduction in detection efficiency was observed. These findings offer invaluable insights into the intricate 
interplay between stellar activity and astrometric precision, significantly advancing our understanding 
in the search for habitable planets. 
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detection 


1. INTRODUCTION 


Over several decades, the discovery of more than 5500 
exoplanets has been achieved through the application 
of diverse methods '. The predominant techniques 
employed for these discoveries are the transit and ra- 
dial velocity (RV) methods. However, these two ap- 
proaches are susceptible to the impact of stellar activ- 
ity, capable of inducing RV variations akin to those at- 
tributed to orbiting planets. With the development of 
the new generation of extreme-precision spectrographs, 
the radial-velocity (RV) measurement at sub-meter-per- 
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second precision becomes feasible, including the Echelle 
SPectrograph for Rocky Exoplanets and Stable Spectro- 
scopic Observations (ESPRESSO) (Pepe et al. 2010), 
the EXtreme PREcision Spectrograph (EXPRES) (Ju- 
rgenson et al. 2016) and NEID (Lin et al. 2022). The 
previous investigations showed that the amplitude of RV 
variations resulting from stellar activity ranges in the 
order of a few m/s (Hatzes 2002; Lagrange et al. 2010; 
Butler et al. 2006, 2017; Laliotis et al. 2023). The disper- 
sion in radial velocity due to granulation is estimated to 
be larger than 0.3 m/s, significantly surpassing the sig- 
nals emanating from terrestrial planets (approximately 
~ 0.09 m/s). Zhao et al. (2022) reported that EXPRES 
is unable to detect Earth-like planets in the presence of 
stellar variability. Periodic variability in stellar activity 
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can mimic a planetary transit signal, leading to potential 
false positives (Lanza et al. 2009). Photometry observa- 
tions are also influenced by stellar oscillations and gran- 
ulation, which are comparable to the planetary signal 
at the precision of CHEOPS (CHaracterising ExOPlanet 
Satellite) and ESPRESSO (Sulis et al. 2023). The simu- 
lations of PLATO (PLAnetary Transits and Oscillations 
of stars) concluded that the noise due to granulation is 
about 100 ppm (Morris et al. 2020), comparable to the 
depth of Earth-like planets. Despite efforts to bolster 
the performance of the RV method through denoising 
techniques such as Gaussian processes and extended ob- 
servational durations (spanning ten years), uncertainties 
in mass estimations persistently exceed 30% (Meunier 
et al. 2023). 

'The astrometric method offers distinct advantages in 
the detection of terrestrial planets. Several astrometric 
missions, among them the Gaia satellite, have made exo- 
planet detection a primary scientific objective. Notably, 
the discovery of over 60 Jupiter-mass planets through 
astrometry was reported in Gaia DR3 (Gaia Collabo- 
ration et al. 2023). Numerous upcoming missions are 
specifically dedicated to the search for Earth-like plan- 
ets, including Theia mission (The Theia Collaboration 
et al. 2017), the Closeby Habitable Exoplanet Survey 
(CHES) mission (Ji et al. 2022), Small-JASMINE mis- 
sion (Kawata et al. 2023), and the Nancy Grace Roman 
Telescope (Gandhi et al. 2023), Habitable Worlds Ob- 
servatory (HWO) survey (Mamajek & Stapelfeldt 2024). 
For instance, Theia endeavors to identify approximately 
six habitable terrestrial planets, a breakthrough with 
profound implications for understanding the conditions 
conducive to planetary formation and the emergence of 
life (The Theia Collaboration et al. 2017). 

'The meticulous consideration of stellar activity noise 
proves indispensable when operating at ultra-high pre- 
cision. Predominantly, the primary astrometric shifts 
arise from the intricate interplay of stellar spots and 
faculae. Stellar spots, magnetic structures that man- 
ifest as darkened regions on the stellar surface, com- 
prise a darker core called the umbra, surrounded by 
a slightly bright halo known as the penumbra (Solanki 
2003). Faculae, existing in proximity to spots, exhibit a 
slightly higher temperature than the photosphere. Con- 
sequently, while facular contrast in temperature is less 
pronounced compared to spots, their significantly larger 
areas contribute crucially to astrometric shifts. 

Considerable investigation has been dedicated to the 
intricate task of modeling astrometric jitter induced by 
spots, particularly in anticipation of the Space Inter- 
ferometry Mission (SIM) (Eriksson & Lindegren 2007; 
Makarov et al. 2009, 2010). Their findings suggested the 


feasibility of detecting Earth-like planets around quies- 
cent stars analogous to the Sun. However, challenges 
arise when dealing with more active stars. Lagrange 
et al. (2011) extended their considerations to encompass 
both plages and spots, utilizing solar spots groups and 
bright structures from USAF/NOAA and MDI/SOHO 
and revealed that the amplitude of activity-induced as- 
trometric signals remains below 0.2 pas. 

Based on their model, Meunier et al. (2019) incorpo- 
rated granulation and chromospheric emission, extend- 
ing the spectral type to F6-K4 stars. They determined 
that the root mean square (rms) of the astrometric ac- 
tivity time series for other main-sequence stars can be 
two to five times the solar value. Nevertheless, the detec- 
tion rates for habitable terrestrial planets persist above 
5096, marking a substantial improvement compared to 
the RV method (Meunier et al. 2020). In their analysis of 
simulated time series, Meunier & Lagrange (2022) con- 
centrated on 55 nearby stars in the Theia sample. Their 
findings echoed similar conclusions, suggesting that stel- 
lar activity minimally impacts the detection of planets 
for solar-type stars, except for the closest stars o Cen 
A and B. They also emphasized the need to enhance 
the observation strategy for certain subgiants due to the 
more distant habitable zone. Additionally, studies based 
on Gaia accuracy, such as Morris et al. (2018), developed 
a model attributing stellar jitter to spots and concluded 
that Gaia’s precision is adequate to measure photocen- 
ter jitter for some M dwarfs within 10 pc. 

In a recent study, Shapiro et al. (2021) computed the 
displacements observable in both the Gaia and Small- 
JASMINE passbands (Kawata et al. 2023). Notably, 
they identified a jitter amplitude in the G-band of ap- 
proximately 0.5 yas for Sun located at 10 pc, antici- 
pated to be twice as large as that in the Small-JASMINE 
near-infrared passbands. Expanding their investigation 
to stars observed at diverse inclinations, metallicity lev- 
els, and active-region nesting degrees (Sowmya et al. 
2021), they discovered that astrometric shifts peak at 
an intermediate inclination in the Gaia-G band. Build- 
ing upon the understanding that active magnetic struc- 
tures tend to emerge in higher latitude regions (Strass- 
meier 2009), Sowmya et al. (2022) advanced their re- 
search by exploring stars with solar fundamental param- 
eters but faster rotation than the Sun. Their findings 
indicated that magnetic jitter is predominantly spot- 
dominated for rapid rotators. In an exploration of the 
multi-wavelength strategy in astrometry, Kaplan-Lipkin 
et al. (2022) employed solar simulations across five dis- 
tinct passbands based on Shapiro et al. (2021). Their 
results highlighted that the minimum detectable mass, 
through a combination of observations like Gaia G and 
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R bands, is 0.005 Earth mass. However, it is crucial 
to note that their assumptions involved perfect instru- 
ments, which can introduce much larger noise than stel- 
lar activity in practical scenarios. 

The primary goal of this work is to evaluate the impact 
of stellar activity on the astrometric signal of the target 
stars of CHES. We developed a model that incorporates 
spots and faculae to simulate stellar activity in astrome- 
try. The stellar activity levels were derived from photo- 
metric data obtained from TESS and the chromospheric 
index logi. Subsequently, we generated simulated 
astrometric time series for each target. Noteworthy, 
we observed that nearly 7596 of stars exhibit minimal 
photocenter jitter, with standard deviations lower than 
0.24 uas. Moreover, we introduced simulated habitable 
planets and conducted retrieval experiments to ascer- 
tain their mass and period. The retrieval performance 
proved to be satisfactory, with approximately 8096 of 
simulated planets around nearly 95% of the stars being 
successfully recovered in our assessments. These find- 
ings underscore the feasibility of the astrometry method 
for detecting terrestrial planets, even when considering 
the noise introduced by stellar activity. 

'The paper is structured as follows: Section 2 provides 
a brief description of the properties of targets and plan- 
ets. In Section 3, we elaborate on the stellar activity 
model, delineate the properties of spots and faculae, and 
elucidate the procedure for generating simulated stellar 
jitter at each time step. The methodology employed 
to estimate the detection efficiency of stars is outlined 
in Section 4. Subsequently, Section 5 presents the out- 
comes of simulations and the detection efficiency. In the 
final section, we summarize the major results. 


2. TARGETS AND APPROACH 
2.1. CHES satellite 


The Closeby Habitable Exoplanet Survey (CHES) 
mission is set forth to unveil habitable-zone Earth-like 
planets surrounding solar-type stars in close proxim- 
ity, approximately 10 parsecs distant from our solar 
system with micro-arcsecond relative astrometry tech- 
niques. The primary scientific objectives of CHES en- 
compass the detection of Earth Twins or terrestrial plan- 
ets within habitable zones orbiting 100 FGK nearby 
stars (Ji et al. 2022). Moreover, the mission aims to ex- 
ecute an exhaustive survey and intricately characterize 
the neighboring planetary systems, and reveal the real 
planetary masses and 3-dimensional orbits that provide 
crucial clues to their planetary formation and dynami- 
cal evolution (Youdin & Goodman 2005; Johansen et al. 
2007; Jin et al. 2016; Huang et al. 2019; Liu et al. 2019; 
Pan et al. 2022, 2024). CHES is expected to enrich the 


sample of Earth-like planets, which helps to comprehen- 
sively understand the formation and evolution of plane- 
tary systems (Ida & Lin 2004; Mordasini et al. 2012; Jin 
et al. 2014; Liu & Ji 2020; Huang & Ji 2022; Huang et al. 
2023). The CHES satellite, outfitted with a 1.2-meter 
aperture high-quality telescope featuring low distortion 
and high stability, will be stationed at the Sun-Earth 
L2 point, continuously observing all target stars over 
a 5-year duration (Ji et al. 2022). Within a planetary 
system, the host star undergoes a subtle wobble around 
the common center of mass, typically demonstrating an 
amplitude of 0.3 pas, which is measured for a solar-mass 
star situated at a distance of 10 parsecs, influenced by an 
Earth-mass planet located at 1 AU. In the astrometry, 
our measurements are focused on the position of the stel- 
lar photocenter, operating under the assumption that it 
coincides with the center of mass. However, beyond the 
gravitational influence of planets, the stellar photocenter 
can experience displacement due to the magnetic activ- 
ity of the host star. Bright faculae regions contribute to 
increased brightness and attract the photocenter, while 
dark spots exert the opposite effect. The resulting pho- 
tocenter jitters are measured in milli stellar radii, a scale 
comparable to the astrometric signature of an Earth-like 
planet (Catanzarite et al. 2008). 


2.2. Target Selection 


In this study, we have identified 94 main sequence 
stars from the CHES targets catalogue. The distribution 
of effective temperature and luminosity is illustrated in 
Figure 1, with colors denoting the stellar types (e.g., 
F, G, K). Remarkably, the majority of these stars are 
positioned within approximately 10 pc from the solar 
system. Additionally, the astrophysical parameters of 
the target stars are derived primarily from the TIC8.2 
input catalogue (Paegert et al. 2021) and are succinctly 
presented in the appendix. 


2.3. Planetary signals 


'The primary scientific goal of CHES is to detect the 
habitable planets orbiting nearby solar-type stars by ob- 
serving the stellar wobbles induced by planetary gravity. 
'The amplitude of the astrometric signal produced by a 
planet can be described as (Perryman et al. 2014): 


mos (m) (n) (e) (rm) e0 


where Mp denotes the planetary mass, M. is the stel- 
lar mass, a represents the semi-major axis of the planet, 
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Figure 1. The Hertzsprung-Russell diagram depicts CHES 
targets, comprising 19 F stars, 24 G stars, and 51 K stars. 


and D is the stellar distance to the solar system. Con- 
sequently, the astrometric signal induced by an Earth- 
mass planet at 1 au, causing a solar-mass star to wobble 
10 pc away, is measured at 0.3 uas as aforementioned. 


3. STELLAR ACTIVITY MODEL 


We formulated a model to simulate the emergence of 
spots and faculae on solar-type stars. Subsequently, we 
generated synthetic observational data, encompassing 
both photometry and astrometry, spanning a five-year 
duration to match the expected time baseline of CHES. 


3.1. Analytic centroid approximation 


The difference in intensity of stellar magnetic features, 
including spots and faculae, compared to the photo- 
sphere, leads to variations in the stellar photocenter as 
the star rotates and features evolve. To evaluate the in- 
fluence of jitters on the detection of terrestrial planets, 
we constructed a model employing butterpy (Claytor 
et al. 2022) to simulate the presence of spots and faculae. 
Subsequently, we computed the astrometric time-series. 

We considered a star with a radius R, and defined 
Cartesian coordinates on the stellar surface to calculate 
variations in the photocenter across the stellar surface 
(Morris et al. 2018). The origin is placed at the cen- 
ter of the star, with the x-axis aligned along the stellar 
equatorial plane. The y-axis indicates the rotation axis, 
and the sign of the z-axis determines whether the feature 
is obstructed by the star. Thus, the position of the ith 
feature can be described as r; = (xi, yi), where r; = [ri]. 

The total flux from the star is characterized by (Morris 
et al. 2018) 


Rx 
= f 2rI(r)dr, (2) 


where r is the distance to origin in unit of stellar radius, 
and I(r) is limb-darkening law. We use a four-parameter 
the non-linear limb-darkening law introduced by Claret 
(2000), 


4 
I(u) u 1—- M a(l- Uu?) H > peri 3 
I() 1 (3) 


0 H< leri, 


where u represents the position on the surface (co- 
sine between the local surface and the line of sight), 
u = vVl—r?, Heri is the drop-off point, or are limb- 
darkening coefficients, which are linked to stellar effec- 
tive temperature and logg (Claret 2018). 

We generated a bulk of spots and faculae in the stel- 
lar surface, and the flux contribution from each feature 
is determined by its position, area, and intensity. The 
intensity of the ith spot can be described as (Meunier 
et al. 2019): 


Fp i = —Csp i Ai * Ui, (4) 


where Csp, is the contrast of spot, A; is the area of spot. 

Here we assume that all features are circular, with a 
radius R;, a position r; and a true area Ap. Features 
near limb will be foreshortened, approximated by an el- 
lipse with a semi-major axis R; and a semi-minor axis 
Riy/1-—r?. Thus the projected area can be expressed 
as: 


The projected area is independent of R;, so we only 
concern on true area and position of spots or faculae. 

In this study, the temperature contrast Csp is utilized. 
We assume that both the photosphere and the spot ex- 
hibit blackbody characteristics. Therefore, Cy, is de- 
scribed through the Planck function f (Meunier et al. 
2019): 


Cup = 1 — Zell) (6) 


Tejp) Loi) 
And the intensity of faculae is similar as spots: 
Fiai AiCts* Mi- (7) 


The astrometric contribution in the 2 direction 18 the 
same for faculae and spots: 


R, Fi 
40 = FLR” 
(8) 


R; DN Cm 
= 0.465 x 10? > DU : 
i (zz) (그) SE SS 


the formula in y direction is similar. 
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3.2. Fundamental stellar parameters 


In our model, the astrometric precision is directly in- 
fluenced by the stellar radius in Equation 8, while the 
planetary signal is connected to the stellar mass in Equa- 
tion 1. In addition to mass and radius, the activity level 
of star determines numbers and distribution of magnetic 
features. In this paper, we adopt the dimensionless pa- 
rameter activity level defined by Aigrain et al. (2015), 
with the assumption that the Sun's level is 1. This pa- 
rameter will be discussed further in section 3.3. 

Subsequently, our objective is to determine temporal 
variations in stellar activity, encompassing both the ro- 
tation period and cycle period. Rotation is a common 
observation in many stars, with typical periods ranging 
from a few days to a few hundred days. Due to the 
planned uneven sampling strategy during CHES obser- 
vations (Ji et al. 2022), the shortest cadence is about 
two or three days, and short-term rotation also impacts 
astrometry. The rotation rate of each target is yet to 
be determined for simulating. Previous studies have de- 
duced stellar rotation periods from lightcurves or spec- 
tral observations (McQuillan et al. 2014; Pan et al. 2020; 
Jin et al. 2023), revealing a correlation between rota- 
tion period and characteristics in large star samples from 
TESS or Kepler (Witzke et al. 2020; Reinhold & Hekker 
2020). 

We have chosen to utilize the empirical law that estab- 
lishes a relationship between the Rossby number R, and 
the estimated turnover timescale 7, with the rotation pe- 
riod. The rotation period can be estimated according to 
Noyes et al. (1984): 


Prot = Ro X Te, (9) 


where the Rossby number is estimated by 108 (Ma- 
majek & Hillenbrand 2008), and the turnover timescale 
is linked to B-V (Suárez Mascareno et al. 2016). Figure 
2 presents the theoretical and observed rotation period 
of 33 stars in our samples. The observed periods closely 
align with theoretical values, falling within the margin of 
error for roughly 8096 of the stars. Despite the more sub- 
stantial disparity between the two stars with the shortest 
periods, our attention is directed towards stars charac- 
terized by moderate periods. This arises from the fact 
that only approximately 1096 of FGK stars manifest ro- 
tation periods of less than five days in the Kepler and 
K2 samples (McQuillan et al. 2014; Reinhold & Hekker 
2020). Consequently, we derived rotation periods for 
stars lacking observed periods with Equation 9. 

Stellar activity cycles are widely hypothesized to arise 
from dynamo processes. The solar cycle spans approxi- 
mately 11 years, and there is evidence suggesting a cor- 
relation between activity cycle periods and rotation pe- 
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Figure 2. Theoretical and observed rotation periods of 33 
stars in our sample, the dashed line denotes the observed 
periods equal to theoretical periods. Errorbars indicates the 
period ranges from uncertainties of logRix. 


riods (Noyes et al. 1984; Brandenburg et al. 1998). In 
our study, we employed a relationship tailored for FGK 
stars to estimate the cycle period (Irving et al. 2023). 
'The cycle period and rotation period details are pre- 
sented in the appendix. All target cycle periods are 
shorter than five years (the time baseline of CHES), en- 
abling us to assume that the maximum stellar activity 
jitter in a cycle is estimated in our simulations. 


3.3. Determination of activity level 


Given that the activity level serves as an adjustable 
scale factor governing the average rate of spot and facula 
emergence (Aigrain et al. 2015), it becomes imperative 
to ascertain the activity level for each target. We de- 
termine the activity level using the chromospheric index 
logR'y and photometric variability as detailed below. 

Initially we derived the lightcurves of 78 stars observed 
by TESS by the package lightkurve (Lightkurve Col- 
laboration et al. 2018), encompassing the simple aper- 
ture photometry (SAP Flux) and the presearch data 
conditioning SAP flux (PDCSAP Flux). For both of 
type data, we removed outlier and bad data based on 
quality flag of TESS. Each light curve was scrutinized 
visually to avoid systematic errors. And the SAP Flux 
is preferred as the long-term trends are removed in PD- 
CSAP Flux. Subsequently, the peak-to-peak amplitude 
values were derived to characterize variability of each 
lightcurve (Hojjatpanah et al. 2020). Given that the 
timescale of a single sector is much shorter than the 
stellar cycle period, for stars observed in two or more 
sectors, we individually processed each light curve and 
calculated the peak-to-peak values. In most cases, the 
differences in peak-to-peak values between different sec- 
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tors are quite small, and the mean peak-to-peak values 
are derived as input parameters in our model. For indi- 
vidual stars, we found that the peak-to-peak value from 
one sector is significantly higher than the others, and no 
systematic errors were found. In this case, the lightcurve 
of this sector is considered to have been observed ob- 
served during the more active phase of the star, and the 
mean value of PDCSAP Flux and SAP Flux of was used. 
Additionally, we calculated the variations of Sun from 
25 Feb. 2003 to 25 Feb. 2020 as a comparison. 

We identified a quadratic relation between peak-to- 
peak values and 10877, which represents the contri- 
bution of the Call H and K lines to the bolometric lumi- 
nosity of the star (Noyes et al. 1984). This index can be 
derived through spectral observations, such as those ob- 
tained by instruments like LAMOST (Zong et al. 2018). 
Figure 3 displays Jog Hun and peak-to-peak variabil- 
ity of 78 samples and Sun (left panel) The variation 
of the Sun was calculated based on Total Solar Irradi- 
ance(TSI) data '. Photometric variations were utilized 
to relate logR'y g activity level in our model. We estab- 
lished a grid of levels from 0.5 to 5, and generated 100 
simulated lightcurves at each level. Figure 4 reveals a 
parabolic relationship between activity level and peak- 
to-peak variation, consistent with the results of Aigrain 
et al. (2015). Consequently, activity levels of targets 
with observed lightcurves could be determined. 

For stars not observed by TESS or Kepler, we fit- 
ted the relationship of activity levels and logR'yg to 
ascertain their activity levels. Figure 3 depicts relation- 
ship between logR'yg and activity level (right panel). 
We found there is a quadratic relation between activity 
level and logR'y g in our samples and the Sun, and the 
relations for stars of different spectral types (FGK) are 
consistent. Notably, the star with the highest activity 
level is HD 17925, an RS CVn Variable star considered 
an outlier and excluded during the fitting process. We 
infer that this relation is applicable to other solar-type 
stars as well. 


3.4. Simulated spots and faculae 


'The primary parameter categories of spots and faculae 
include distribution, variability, contrast and size. We 
utilize the Python package butterpy to generate spots 
with spatial and temporal distribution (Claytor et al. 
2022). 

Over the approximately 11-year solar activity cycle, 
spots emerge within active regions at high latitudes ini- 
tially. Subsequently, these regions migrates toward the 


1 TSI data is from observation of The Solar Radiation and Climate 


Experiment (SORCE) 


equator, indicating some overlap between consecutive 
cycles, resulting in a butterfly diagram illustrating the 
latitude distribution of spots over time. In our model, 
we set the overlap duration to be 10% of cycle period. 
The butterfly diagram has also been observed in other 
stars (Bazot et al. 2018; Nielsen et al. 2019). In our 
model, we apply the butterfly spot emergence feature of 
butterpy. Because the faculae typically existed around 
spots, we add a facula near each spot in our model 
(Borgniet et al. 2015). The input parameters needed for 
butterpy are given in Table 1. Active regions are posi- 
tioned at medium latitudes in our samples. Since polar 
spots have been observed in several fast rotators using 
Doppler and Zeeman Doppler imaging (Vogt & Penrod 
1983; Strassmeier 2009), we also generated spots and 
faculae at latitudes ranging from 55° to 85° for samples 
with periods shorter than five days. We found that the 
difference in standard deviations of astrometric jitters 
between these two models (i.e., medium latitudes and 
high latitudes) is approximately 0.1 yas. The latitudes 
of the activity regions are considered less influential due 
to the symmetry of our model. 

We estimated the temperature difference between the 
photosphere T, and spots Tsp of FGK stars (Berdyug- 
ina 2005; Namekata et al. 2019), revealing a trend with 
effective temperature T; of the star: 


ANTT, 
=3.58 x 10-5 T, + 0.249T, — 808 (K). (10) 


As the area ratio between penumbral and umbral 
of sunspots is approximately 4 to 5 (Solanki 2003), 
the penumbral is significantly hotter than the umbra 
(~ 1000K) (Johnson et al. 2021). Therefore, we weigh 
the differential temperature based on the relative areas 
in the penumbrae and umbrae within the range 0.2 and 
0.6. The coefficient used is similar to that employed by 
Lagrange et al. (2010). The contrasts of stellar faculae 
are poorly constrained, we use an intensity contrast in 
the optical band for G-type stars (Borgniet et al. 2015): 


Cia = 0.131618 — 0.218744, + 0.1047572, (11) 


the simulated spectra show that the differences in facu- 
lar contrasts between K-type and G-type stars are rel- 
atively small at the band of CHES, and the contrasts 
of F-type stars are slightly lower than GK stars (Nor- 
ris 2019; Norris et al. 2023). The contrasts of faculae 
are much lower than those of spots and are observed 
generated near spots. However their area is greater as 
observed on the Sun, the contribution from faculae is 
therefore expected to be significant. In this paper we 
use the area ratio to describe the size of faculae because 
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Figure 3. Left panel: the relation between logR'yg and peak-to-peak variations of photometry, the dashed line shows a 
parabolic fitting. Right panel: same for the relation between log y and activity level, while errorbars shows the uncertainties 
of activity level. The blue circle, orange cross, and green triangle correspond to F,G,K stars respectively. 
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Figure 4. The activity level with photometric peak-to-peak 
variations. The step of activity level is 0.2, the y-coordinate 
of each point indicates the mean peak-to-peak variations of 
100 simulated lightcurves. 


we add a facula near each spot (Borgniet et al. 2015), the 
overlap situation was not considered in our simulations, 
given that the areas of features are significantly smaller 
compared to the stellar disk. We employed a lognor- 
mal distribution for spots in our simulations (Bogdan 
et al. 1988; Solanki & Unruh 2004; Baumann & Solanki 
2005). Additionally, we used a lognormal distribution to 
describe the ratio between faculae and spots (Borgniet 
et al. 2015). 

Given that the emergence of spots and facula occurs 
within a few hours (Howard 1992), smaller than our time 
step, and we focus here on planets with periods of a 
few hundred days, we neglected the growing process and 
assumed that all features are created with maximum 
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Figure 5. An example of the temporal and spatial evolution 
of spots and faculae on the surface of target 107 Psc, with a 
cycle period about 1200 days. Stellar magnetic features lat- 
itude distribution follows the butterfly diagram during eight 
years, the circles represent spot, and the crosses represent 
facula. The colorbar shows logarithmic areas of features, in 
unit of MSH (millionths of solar hemisphere). 


area (Borgniet et al. 2015). We employed a parabolic 
decay law for both spots and faculae (Gilbertson et al. 
2020). 

All input parameters of spots and faculae for our 
model are summarized in Table 1. For each star, the 
simulations last for eight years with a 0.9 day time step 
(total 3000 points) to ensure the simulation covers the 
entire cycle. At the beginning of each simulation no 
features are present on the stellar surface. Thus, we ex- 
tracted the middle five years, assuming the simulation 
has reached a steady state. Figure 5 displays an exam- 
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Table 1. Input parameters for spots and faculae 


Category Parameter Value/range Unit Reference 
Stellar activity Cycle Period 0-60 yr 
Cycle Amplitude 0.1-5 - 
Cycle Overlap 0.1 - Aigrain et al. (2015) 
Spots properties Mean initial size 46.5 MSH” 
Standard size deviation 2.14 MSH 
Max size 1500 MSH Meunier et al. (2019) 
Decay rate 43.9 MSH/d 
Max ave latitude 35 deg 
Min ave latitude 7 deg Claytor et al. (2022) 
Faculae properties log Mean size ratio 0.8 - 
log Standard size ratio 0.4 - 
Latitude deviation 3 deg 
Longitude deviation 3 deg 
Decay rate 21.9 MSH/d  Borgniet et al. (2015) 
NOTE—a:micro stellar hemispheres 
ple pattern of 107 Psc to show stellar magnetic features 
latitude distribution over eight years. 
0407 . @ Spots 
x  Faculae 
4. DETECTION EFFICIENCIES 0.08 - 
Moreover, we used the inject-recover method to esti- doel 
mate detection limits of each star. We randomly gen- " e 
erated 5000 planets in the parameters space. Since 2 Wir x X M 
the main goal of CHES is to search Earth-like planets, ub e ^ XY e" R 
we set the planetary mass range to of 0.5 to 5 Earth ra eg? vi "E 
masses. All planets were placed in the habitable zone, 0.00 F n ^: ox ti 
with the inner and outer semi-major axis boundaries BH toe Sey T 
calculated based on the Runaway Greenhouse limit and e x o is 
Maximum Greenhouse limit, respectively (Kopparapu -oo4[ | ! | | $9 Fe. 
et al. 2014). All orbital elements and planetary masses x M ME E E t n 0:92. 50.04 — .0:06 


followed a uniform distribution. The stellar wobbles due 
to a planet represent Keplerian orbit motion and can 
be derived through the Thiele-Innes equations (Thiele 
1883). Apart from the photocenter shifts due to stellar 
activity, we also account for other measurement errors, 
including photon noise, telescope error, and calibration 
errors. Gaussian noise, with a mean of 0 and standard 
deviations of 0.7 pas is added in the simulations (Ji et al. 
2022) . 

For planets retrieval, the Lomb-Scargle periodograms 
are utilized (Horne & Baliunas 1986). The period with 
the highest power is considered as a prior during fit- 
ting. Then we use a sinusoidal function based on the 
Thiele-Innes equations to fit planetary signal (Jin et al. 
2022). In the fitting process, we assume all planets are 
in circular orbits to reduce computation. The ampli- 
tude of the sinusoidal function is denoted as o, and the 
planetary mass is derived using Equation 1. We con- 


Figure 6. The photocenter jitter due to spots and faculae 
of HD 103095 during five years, the blue circles and orange 
crosses represent effects of spots and faculae, respectively. 


sider a planet to be recovered if both the derived period 
and mass within 3096 error is recovered. Additionally, 
we use the SNR criteria as a comparison. The SNR in 
astrometric planets is defined as (Unwin et al. 2008): 


(12) 


c 
where the c is the standard deviations of total noise, 
Nyis is the number of observation, depending on obser- 
vational strategy of CHES (Tan et al., in prep.). Table 
2 further presents DHZ (the distance to the center of 
habitable zone) for each target. 
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5. RESULTS OF JITTER AMPLITUDE AND 
DETECTION EFFICIENCY 


We applied the methods outlined in Section 3 to 94 
stars. First we found there is a quadratic relation be- 
tween activity level and the photometric peak-to-peak 
variations, as expected by Aigrain et al. (2015). We 
then fitted a parabolic relation to calculate activity lev- 
els from og Ru for our samples, and as the Sun corre- 
sponds to this relation as well, we infer that the relation 
can be extended to other FGK stars. Subsequently, we 
generated astrometric timeseries due to spots and fac- 
ulae for each target. As an example, we consider the 
target HD 103095, a K-type star with no detected plan- 
ets. The astrometric jitters of its photocenter are shown 
in Figure 6, with an amplitude about 0.1 uas, which is 
lower than signal from an Earth-mass planet (~ 0.3 pas). 
We derived the standard derivations in two directions of 
all targets, as shown in Figure 7. It was found that for 
more than 95% of the targets, the standard deviations 
of photocenter shifts due to activity are less than 1 pas. 
Figure 8 shows the mean shifts of photocenter in our 
samples, which are within 0.4 pas for all stars. And the 
errorbars represents amplitudes of jitters. We found sev- 
eral stars with amplitudes larger than 1 uas, because of 
shorter distance (e.g. Alpha Cen A and eps Eri) or the 
combination of higher activity level, larger size and short 
distance (e.g. eta Cas A). We found the jitter in the Y 
direction is slightly smaller than that in the X direction 
because of the symmetry of butterfly diagram. And we 
derived the detection rates corresponding to terrestrial 
planets in habitable zone. The results for all targets are 
provided in Table 2. The inject-recover results of 107 
Psc is shown in Figure 9. We found about 95% planets 
with SNR > 6 were recovered, consistent with those of 
Perryman et al. (2014). 

'The histogram of detection efficiency is shown in the 
left panel of Figure 10, with nearly 9596 stars having a 
detection rate above 8096. The detection efficiency de- 
creases for some cold stars, such as HD 32450 and HD 
21531, due to the weak planetary signal from a closer 
habitable zone. Despite several close stars, such as « 
Cen A and eps Eri having an obvious jitter, their detec- 
tion efficiencies remain high because the planetary sig- 
nals also increase as the distance gets closer. The right 
panel of Figure 10 displays the detection efficiencies ver- 
sus the distance of the stars. In our sample, six stars 
have detection efficiencies below 8096; five of them are 
located 8 pc away, exhibiting a slightly weak planetary 
signal. Notably, the subgiant alf CMi has a detection 
efficiency of about 79.896 due to its much more distant 
habitable zone. The period corresponding to habitable 
planets around alf CMi is approximately the same length 
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Figure 7. The standard deviations of astrometric jitter in 
two directions of all targets from simulations. The blue circle, 
orange cross, and green triangle correspond to F,G,K stars 
respectively. The two histograms shows jitter’s distribution 
in our samples. 
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Figure 8. The mean jitter caused by stellar spots and fac- 
ulae from simulations. The errorbars indicate the maximum 
and minimum jitter duration five years. 


as the time baseline of CHES, making planet character- 
ization difficult. 


6. CONCLUSIONS AND DISCUSSIONS 


In this study, we conducted simulations to evaluate 
the influence of stellar activity on the detection of ter- 
restrial planets in the context of the ultra-high precision 
astrometry space mission, CHES, which aims to perform 
measurements on approximately 100 nearby stars. Our 
model, implemented within butterpy, generated simu- 
lated astrometry data for the identified CHES targets. 
We estimated the contrast of spots based on Equation 
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Figure 10. Left: The histogram of detection efficiency, the colorbar show mean effective temperature in each bin. Right: The 
detection efficiency and distance of all targets, sizes of each point indicate stellar radius. The blue circle, orange cross, and green 


triangle correspond to F,G,K stars respectively. 


10, which applies to stars with effective temperatures 
between 3000 to 7000 K (Berdyugina 2005; Namekata 
et al. 2019). Notably, there are about a dozen stars be- 
yond this temperature range. As the empirical formula 
is believed to overestimate contrasts, especially for hot- 
ter stars (Borgniet et al. 2015; Meunier et al. 2019). 
Borgniet et al. (2015) also estimated a lower bound of 
spots contrasts. In this paper, we choose Equation 10 
as an extreme case because higher contrasts cause larger 
astrometric noise. Similarly, the facular effect of F-type 
stars can also be overestimated in our simulation as in- 


troduced in Section 3.4. Therefore, higher detection ef- 
ficiency of CHES is actually possible. 

As some peak-to-peak values were derived from multi- 
ple sectors, there are errors in each star’s activity level, 
as depicted in Figure 3. The astrometric jitters of ten 
targets with larger errorbars were calculated based on 
their maximum activity level value. Tet Per exhibits a 
noise about 1.3 pas with an activity level of approxi- 
mately 2.5. While noise from stellar activity is still not 
the dominant factor for other stars, with standard devi- 
ations lower than 0.7 pas. These ten stars are indicated 
in Table 2. 


ASTROMETRIC NOISE AND PLANETARY DETECTION EFFICIENCY DUE TO STELLAR SPOTS AND FACULAE 1l 


'The determination of each star's activity level relied 
on the chromospheric index log; g. Our observations 
revealed a discernible correlation between increased ac- 
tivity levels and heightened variability in both photom- 
etry and astrometry, best described by a parabolic rela- 
tion. Subsequent analysis of the time series facilitated 
the estimation of detection efficiencies for planets within 
the habitable zone of each target. Notably, the influence 
of stellar activity on astrometry was found to be mini- 
mal, with photocenter jitter measuring less than 1 pas 
for around 9596 of the targets. In the context of terres- 
trial planet detection, the astrometric method demon- 
strated superior efficacy compared to current radial ve- 
locity and transit methods. 

We compared our findings with the results of Meunier 
& Lagrange (2022), who estimated the astrometric noise 
from stellar activity of Theia's targets. Among the 37 
common targets, the differences in astrometric jitter are 
less than 0.3 pas for approximately 80% of the stars. 

In this work, we intentionally exclude the considera- 
tion of other stellar characteristics, such as oscillation, 
granulation, and supergranulation, as their expected im- 
pact on astrometry is deemed negligible (Sowmya et al. 
2021). The associated timescales (Staude et al. 1999; 
Borgniet et al. 2015; Chiavassa et al. 2017) are consider- 
ably shorter than the observation cadence of the Closeby 
Habitable Exoplanet Survey (Ji et al. 2022). Our focus 
remains on single-planet systems, although multi-planet 
systems are prevalent in observational datasets. Ap- 
proximately 10 targets in our sample are identified as 
hosting hot Jupiters or Neptunes based on radial veloc- 
ity measurements. The presence of existing gas-giants 
may obscure signals from Earth-like planets, a factor we 
plan to address in forthcoming study. 

In our simulations, stellar inclinations are consistently 
set in an edge-on configuration. While Meunier et al. 
(2019) concluded that the amplitude of radial velocity 
reaches a maximum in this configuration, Sowmya et al. 
(2021) found that the largest astrometric jitter occurs at 
intermediate inclinations. A nuanced exploration of the 
impact of inclinations is scheduled for future investiga- 
tions. In this work, we do not consider the growth phase 


of activity regions, which is not well understood observa- 
tionally at present. Since the growth rates of sunspots 
are much faster than their decay rates (Howard 1992; 
Forgács-Dajka et al. 2021), and nearly 9096 of the tar- 
gets in our sample exhibit higher activity levels than 
the Sun, we chose to disregard this effect in our study. 
The total flux contribution, including the growth pro- 
cess, is estimated to be less than 1.096 based on simu- 
lations by Aigrain et al. (2015). Assuming that the dip 
in flux originates from only one spot on the edge, the 
photocenter shifts are approximately 0.4 jas for a solar- 
like star at a distance of 10 pc, which is comparable to 
the results obtained from our no-growth model. Given 
the limited understanding of spot evolution, this growth 
phase will be incorporated in future study to provide 8 
more comprehensive description of the astrometric noise 
from stellar activity at small timescales. This approach 
aims to contribute to a more holistic understanding of 
the challenges and opportunities associated with astro- 
metric detection methods, especially when considering 
diverse planetary systems. Furthermore, our forthcom- 
ing research will systematically investigate additional 
factors influencing the astrometric precision of CHES, 
including comprehensive considerations such as the po- 
tential presence of planets around the reference stars and 
the stability of the satellite. 
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